Pattern Stability and Trijunction Motion in Eutectic Solidification 
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We demonstrate by both experiments and phase-field simulations that lamellar eutectic growth can 
be stable for a wide range of spacings below the point of minimum undercooling at low velocity, 
contrary to what is predicted by existing stability analyses. This overstabilization can be explained 
by relaxing Cahn's assumption that lamellae grow locally normal to the eutectic interface. 



The solidification of eutectic alloys is both a striking 
example of spontaneous pattern formation in nature and 
a metallurgical problem of widely recognized practical 
importance jjj. This growth process has been tradi- 
tionally studied by directional solidification experiments 
where a sample containing a binary alloy of near-eutectic 
composition is pulled with a fixed speed V in an ex- 
ternally imposed temperature gradient G. This setup 
produces a wide range of microstructures of which the 
simplest is an array of lamellae of two coexisting a and /3 
solid phases growing into the metastable liquid, as shown 
in Figure |l|a. The steady-state growth of a perfectly pe- 
riodic lamellar array is described by the classic Jackson- 
Hunt (JH) theory H that predicts the relationship 




AT(A) =T E - Tau (A) = K Y VX + K 2 /\, 



(1) 



between the lamellar spacing A (width of one lamella 
pair) and the undercooling, AT(A), which is the differ- 
ence between the eutectic temperature, Te, at which the 
three phases (a, (3 and liquid) coexist in equilibrium, 
and the average temperature of the nonequilibrium eutec- 
tic interface spatially averaged over one spacing, T av (X). 
The first and second terms on the right-hand-side of Eq. 
|l| represent the undercooling necessary to drive the dif- 
fusive transport of the two chemical components of the 
alloy in the liquid (coupled growth) and the capillary 
undercooling associated with the curvature of the solid- 
liquid interface, respectively; K\ and K 2 are constants 
that only depend on the alloy system and the overall 
composition of the sample. 

The JH result implies that the growth undercooling has 
a minimum value AT m — (AK\K 2 ) X I 2 P 1 / 2 for a spacing 



FIG. 1. Photographs: Lamellar-eutectic fronts of a nearly 
eutectic CBr4-C2Cl6 alloy in directional solidification (the 
growth direction is upward) in 1 2-pan thick samples. Graphs: 
Interlamellar spacing A (thin lines) and position z of the front 
(thick lines) as functions of the space variable x. a) Station- 
ary pattern (V = 0.5 fims' 1 ; G = 80K cm' 1 ), b) Modulated 
pattern (V = 0.25 ^ms' 1 ; G = ASKcm' 1 ). 
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known to limit the array stability at large spacing 
JH have credited Cahn in Ref. Q for pointing out 
A c should be equal to A m if one assumes that lamel- 
;row locally normal to the envelope of the eutectic 
face. Langer later formalized this result by showing 
a large-scale and small-amplitude spacing modula- 
of a steady-state array obeys the diffusion equation 



d t X(x,t) =Ddl\(x,t), 



(3) 



X m = {K 2 /K 1 ) 1 '' 2 V 



-1/2 



(2) 



where x is the coordinate perpendicular to the growth 
axis z, and D = D± with 



which is easily found by setting dAT/dX = 0. 

Lamellar growth is well known to be unstable for spac- 
ings smaller than a critical value A c . This instability leads 
to the local elimination of lamellae and is the mechanism 
by which the array increases its average spacing during 
the dynamical transient that produces the final pattern. 
Hence, it is crucially important for understanding pat- 
tern selection in this system. Oscillatory instabilities are 
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We have defined A = A /A m where A is the spacing of 
the steady-state array being perturbed, and we have used 
the subscript "_L" to stress that this expression for D is 
obtained under Cahn's assumption that lamellae grow 
normal to the interface. Langer's analysis reproduces 
the JH-Cahn result that growth is unstable below A m 
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since perturbations are amplified (decay) when D± < 
(> 0). In addition, it shows that lamella elimination is 
initiated by a long-wavelength diffusive instability that 
is generically present in one-dimensional pattern forming 
systems with translation symmetry along x. 

In this letter, we study the steady-state and stabil- 
ity properties of lamellar cutectic growth by thin-sample 
directional solidification experiments in the transparent 
organic system CBr^C^Clg and by two-dimensional sim- 
ulations of a phase-field model, and we extract from both 
approaches independent accurate determinations of A m 
and A c . An important and novel component of our exper- 
iments is the direct measurement of AT(A), which allows 
us to obtain A m from the minimum of this curve rather 
than computing its value from Eq. |^, thereby circum- 
venting uncertainties in materials parameters. We find 
that, in both experiments and simulations, A c is sub- 
stantially smaller than A m , even for typical directional 
solidification growth conditions where the two spacings 
have previously been assumed equal. Furthermore, by 
analyzing the decay of long-wavelength perturbations of 
the array in both experiments and simulations, we obtain 
a direct measurement of £>, which allows us to shed light 
on the origin of the discrepancy between A c and A m . 

The experiments were made with a nearly eutectic 
CBr4-C2Cl6 alloy prepared with zone refined materi- 
als in thin (12 /im thick) glass wall samples (8 mm 
wide and 60 mm long). The values of G used ranged 
from 40 Kcm~ x to 110 Kcm^ 1 (±10%), and those of 
V from 0.125 to 0.75 /xtos _1 (±4%). Details concerning 
the preparation of the samples, the solidification setup, 
and the visualization of the front shape can be found in 
Refs. [@§. 

The steady-state AT(A) curve has never been mea- 
sured directly due to the fact that AT m is usually of the 
order 0.01 K, whereas the absolute temperature is not 
known with a precision better than about 0.1 K. To 
overcome this difficulty, we exploit two key ingredients. 
Firstly, as will be described elsewhere, we are able to cre- 
ate a large-scale modulation of spacing where A(x) varies 
between two extreme values that comprise X m , as shown 
in Fig. |l|(b). Secondly, we measure the z coordinate of 
the solid-liquid interface averaged over one A, which we 
denote by z(x), and compute the local front undercooling 
using T av (x) = Gz(x) + Tq, where To is an unknown con- 
stant. By eliminating x between To — T av {x) and X(x), we 
obtain To — T „(A) which we then fit to Eq. [l] expressed 
in the form T - T av {\) = AT m (A/A m + A m /A)/2 - AT , 
using A m , AT m and ATo = Te — T) as adjustable pa- 
rameters. A plot of To — T av (\) and its fit is shown in 
Fig. H The fit is very good for A smaller than about 
1.25A m . The departure observed beyond this limit is 
compatible with the one which exists between the nu- 
merically calculated curves AT(A) and the JH approx- 
imation ^,0]. We performed such measurements for V 
ranging from 0.125 to 0.5 fims^ 1 . We found X^V — 
K 2 /K 1 = 193 ± 16 /im^- 1 and AT 2 JV = AK X K 2 = 
(2.7 ± 1.3) x 10~ 3 K 2 s/j,m~ 1 . These values compare well 
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FIG. 2. Undercooling To — T av vs spacing A measured ex- 
perimentally from the photograph shown in Fig. lb (open 
circles and dots). Thick line: Best fit of the JH law, Eq. 
to the data represented by open circles. Vertical bar: error 
range. A remarkable feature is the presence of spacings con- 
siderably smaller than A m ~ 27/xm. The value of the smallest 
stable spacing X c predicted from Eq. ^ below is 19 /xm. 

to those calculated from the material constants of CBr-i- 
C 2 C1 6 given in ref. §], namely, \ 2 m V = 185 ±20 /Ws" 1 
and ATl/V = (1.2 ± 0.2) x IQ^K^snrrT 1 . 

To study experimentally the small spacing stability 
limit, we exploit the fact that A m ~ V -1 / 2 . Therefore, 
we can effectively vary Ao/A m by performing downward 
velocity jumps of relatively large amplitude. Namely, 
we start from a stable quasi-stationary periodic array of 
spacing Ao at a higher velocity, and then observe whether 
the same array at the lower velocity, and hence smaller 
Ao/A m , remains stable or becomes unstable. 

To measure experimentally the array diffusion con- 
stant D, we use the fact that the amplitude of a long- 
wavelength modulation of spacing of the form 

X(x, t) w Ao + 5\q e~xjp(ikx + u>f.t) (5) 

decays exponentially in time when the array is stable 
(Ao > A c ). Substituting this form in Eq. ||, we obtain the 
simple dispersion relation u>k = —Dk 2 which is valid if 
the wavelength of the perturbation 2-7r/fc 3> Ao- Knowing 
k, and calculating ujk by a fit of the measured amplitude 
of the modulation to a decaying exponential, we obtain 
D. This is illustrated in Fig. |3| for a case where 2-7r/fc w 
7Ao- We deduce from this measurement that 

D = D ± + D n (6) 

where D» is a positive contribution responsible for the 
overstabilization of the array. The latter is calculated by 
taking the difference between D, and D± evaluated via 
Eq. [| using the value K\ = 1.9 x 10~ 3 Ks/im~ 2 obtained 
from our present experiments. 

Next, we simulate a phase-field model of a two- 
component (AB) eutectic alloy which is completely sym- 
metric under the exchange of a and /3. Our goal here 
is not to model quantitatively the experiments but to 
demonstrate the generality of our results in different 
alloy systems. An order parameter 0, which distin- 
guishes between solid (0 = +1) and liquid (0 = — 1), 
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FIG. 3. Experimental measurements of the spacing A vs 
the space variable x at different times t showing the relaxation 
of a large-scale modulation of a lamellar pattern of Ao ~ 1. 
(G = 80Kcm~ 1 ;V = 0.5 ^ms" 1 ). Inset: Amplitude A of the 
dominant mode (wavelength 145/im) as a function of t, fitted 
by an exponential law (time constant 410s). 



is coupled to the dimensionless concentration field, u = 
(c— ce)/(Ac/2), where c is the mole fraction of B, and Ac 
is the difference between the solid composition in /3 and 
a at Te- The non-conserved and conserved dynamics for 
<j) and u, Td t 4> = -SF/5cj>, and d t u = V • (M{<f>)V5F/5u), 
respectively, are derived from the functional 
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dV 
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|Vu| 
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-|V0| 2 + /(<^,T) 



(7) 



where V is the volume of the system and 

h(<f>) £ , l-h(<P) 



/(</>, u, T) = f D w(4>) 



1 
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fh (8) 



is the bulk free-energy density that is the sum of the 
standard double-well Jdw = —<fi 2 /2 + ( / )4 /4, and a 
concentration-dependent part that interpolates between 
the bulk free-energy of the liquid, /; = u 2 /2, and the 
solid f s = (u 2 - l) 2 /8 - (T E - T)/T E , with h(<j>) = 
3(0 — 3 /3)/2. Furthermore, we used the mobility 
M ((f)) =D;[1-(1 + 0) 4 /16] , where D t is the solute diffu- 
sivity in the liquid, that makes the diffusivity in the solid 
vanish and yields efficient computations 0. Directional 
growth is implemented using the frozen-temperature ap- 
proximation T(z, t) = TE + Gz — Vt, and periodic bound- 
ary conditions in x are used in all simulations. 

Steady-state AT(A) curves were obtained from short 
simulations with two lamellae. The stability was studied 
with long simulations where steady-state arrays of up to 
20 lamellae, constructed from the two-lamella solutions, 
are slightly perturbed by a long-wavelength modulation 
of spacing of the form of Eq. |5|. By exponential fits of 
the amplitude of modulation vs time for different fc, we 
obtain the stability spectrum tok as illustrated in Fig. |l|, 
and hence D by a quadratic fit of to vs fc at small k. 
The stability limit A c is then obtained by determining 
where D changes sign and values of Du are obtained by 
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FIG. 4. Growth rates ujk of spacing perturbations versus 
wave vector k as extracted from phase-field simulations. We 
took W u = Wj, = D t = r = T E = 1 and V = 0.01, which gave 
Xm/lD = 0.333, where Id = Di/V is the diffusion length. We 
measure the strength of the temperature gradient by the ratio 
It /Id, where It = mAc/G is the thermal length (mAc = 2 
in the phase-field units). For both arrays, Ao = 0.84; circles: 
It /Id = 20; squares: It /Id = 2. Dashed lines: fits Uk = Dk 2 
obtained from the points with smallest k. 
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FIG. 5. Dii/(VAo) versus Ao for simulations and experi- 
ments. Dashed line: Eq. bl drawn with A = 0.15. 



subtracting D± from D. We find that the dimensionless 
ratio £)||/(VAo) varies with Ao, but negligibly with G and 
V. Moreover, we find that the simple form 



D||/(VA ) « AAo, 



(9) 



with A « 0.15 gives a reasonable fit to our phase-field 
simulation results as shown in Fig. |^. Remarkably, our 
experimentally determined values of Dm are close to those 
obtained in the phase-field simulations, even though the 
two alloy systems are different. 

To interpret our findings, let us briefly review Langer's 
analysis that yields D = D±. Its first ingredient is the 
assumption that the interface adjusts adiabatically its 
average temperature to the local spacing, or 



AT[X(x,t)] «-GC(M), 



(10) 



where AT(A) is the same as in steady-state and £(x) is 
the z coordinate of the envelope of the eutectic interface, 
defined as a smooth curve interpolating the positions of 
the three-phase junctions (trij unctions), with the origin 
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at Te- The second is Cahn's assumption that lamellae 
grow normal to this envelope, which for a small pertur- 
bation is equivalent to 



d t y(x,t) w -Vd x C(x,t), 



(11) 



where y(x,t) is the lateral displacement of trijunctions 
from their steady-state positions. Finally, it follows from 
the definition of y that X(x,t) sa Ao(l + d x y). Differen- 
tiating both sides of this identity with respect to time 
and using Eqs. |l|, |h| and 10} one obtains the diffusion 
equation (g) with D = D\. 

We checked that Eq. |Hj is indeed faithfully obeyed in 
the range of wavelengths that we consider here. Conse- 
quently, the discrepancy between D and D± must origi- 
nate from a correction to Cahn's normal growth assump- 
tion. It is simple to show that the modified phcnomcno- 
logical form 



d t y(x, t) « -Vd x ((x, t) + DAX(x, t)/X , 



(12) 



yields the diffusion equation (|^) with D given by Eq. |^. 
Equation ^ implies that trijunctions also move locally 
"parallel" to the envelope of the eutectic interface in re- 
sponse to a gradient of spacing. To see physically why 
this lateral motion overstabilizes the pattern, consider a 
local depression in an array of initial spacing Ao- Cahn's 
normal growth assumption implies that such a depres- 
sion will produce a local decrease of spacing, and hence 
a local increase in undercooling that will amplify this de- 
pression if Ao < A m (because dAT/dX < in this case). 
This well-known argument yields A c = A m . In contrast, 
the second term on the right-hand-side of Eq. [l2] implies 
that the lateral motion of trijunctions opposes the local 
decrease in A, and hence helps flatten the interface. 

A prediction for A c = X c /X m can be obtained by set- 
ting D = D± + D|| = 0, which, using Eqs. [|, f| and |, 
yields the cubic equation 



'"if 



AG 
KiV 



A, = 0. 



(13) 



With Ki = 1.9 x W- 3 Ksfim- 2 and A = 0.15, we ob- 
tain A c = 0.70 for the experiment of Fig. [| the lowest 
observed spacings are just above the predicted stability 
threshold, as it should be. For the phase-field simulations 
(A'i = 0.1120), we find A c = 0.942 for l T /l D = 20 and 
A c = 0.715 for It/Id = 2. 

Two previous stability analyses have predicted that A c 
should be smaller than A m . The one by Caroli et al. 
fiof , however, is restricted to a large G limit that cannot 
be compared with our results. The other by Chen and 
Davis |llj does not have this restriction, but predicts a 
departure of A c from X m that is about one order of mag- 
nitude smaller than found here and predicted by Eq. |l3| . 
We believe that the lateral motion of the trijunctions 
is due to a coupling between the diffusion field and the 
non-planar front geometry on the scale of the individual 
lamellae; such effects would appear only at higher orders 



in the analyses cited above. Therefore, an analytical un- 
derstanding of eutectic stability at small spacings from 
sharp-interface models largely remains to be developed. 

The present results shed light on a number of previous 
observations. In metallic eutectics, K\ is generally close 
to 10~ 2 Ksfim- 2 plJl^l, so that the value of V below 
which the departure of A c from unity becomes significant 
is of about 1 /^ms _1 for G in the lOOAcrn -1 range. This 
may explain the deviation from the law A oc U~ ' 5 , where 
A is some empirically defined average eutectic spacing, 
that has been observed at V lower than about 1 /ims -1 
in a number of metallic eutectics PH . Similarly, the over- 
stability due to the lateral motion of the trijunctions may 
explain why coupled growth in a peritectic system has 
recently been found to be stable [g_5| in a situation, anal- 
ogous to that of eutectics at A < A m , where the interface 
should be unstable according to the JH-Cahn stability 
arguments |l6| ] . Finally, our results also improve our un- 
derstanding of the morphological instability that leads 
to the formation of eutectic colonies in the presence of a 
dilute ternary impurity 
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